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Abstract. Consider the model of bipartite entanglement for a random pure 
state emerging in quantum information and quantum chaos, corresponding to 
the fixed trace Laguerre unitary ensemble (LUE) in Random Matrix Theory. 
We focus on correlation functions of Schmidt eigenvalues for the model and 
prove universal limits of the correlation functions in the bulk and also at the 
soft and hard edges of the spectrum, as these for the LUE. Further we consider 
the bounded trace LUE and obtain the same universal limits. 



1. Introduction and main results 

Quantum entanglement has recently been studied extensively [23l [181 1321 123 HSl 
1211 116[ [3] due to its central role in quantum information and quantum computation, 
which is treated as an indispensable resource |20| . The entanglement of random 
pure quantum states is of much interest in the context of bipartite entanglement, 
and statistical properties of such random states are relevant to quantum chaotic 
systems, see [Ml US] and references therein. 

In the present paper, we consider a bipartite quantum system (a system with 
its surrounding environment). Given a composite system A (S) B of an (NM)- 
dimensional Hilbert space ^(^a^) = H^^^ ® n^g\ let |ef)^i and lef)^^^ be 
two complete orthogonal basis states for the subsystems A and B, respectively. 
Without loss of generality, we assume N < M. Any quantum state |$) € A® B 
can be expanded as a linear combination 

N M 

(1-1) |$)=^5]X,,|ef)®|ef) 

i=i j=i 

where these coefficients e C form a rectangular N x M complex matrix X = 
[Xij]. The composite state |$) is fully unentangled (separable) if |$) can be written 
as a direct product of two states l^'^) € A and |$^) e B, i.e., |$) = |$^) (g) 
$'^), otherwise referred it as an entangled state [5D]. The composite state |$) is 
normalized pure state system if the density matrix of 1$) is given by p — |$)($| 
satisfying tr [p] = 1 . The reduced density matrix of the subsystem A by tracing 
over the states of the subsystem B is defined jH Hn] by 

M N 

(1.2) pA = trB[p]=Y.('^f\p\ef}= W.,\et){ef\, 

where Wij are the entries oi N x N square matrix W = XX"^ with tr = 1 due 
to the normalized restriction that tr [p] = 1. Analogously, pb = tr^[/9]. It is not 
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difficult to prove tliat tlie reduced density matrices pA and ps have the same set 
of non-negative eigenvalues xi, X2, ■ ■ ■ , xpf (we call them Schmidt eigenvalues of the 
quantum state |<i>)) satisfying J^iLi^i = 1- Let vf denote the eigenvector of the 
square matrix W corresponding to the eigenvalue Xi. Then p^ can be expressed as 
Pa = J2iLi ^i\''^f){'''f\- A similar representation holds for ps- The composite state 
I $) has a well-known Schmidt spectral decomposition |3] 

N 

(1.3) |<i>)=^^|^;,f)®|^;f). 

i=l 

A pure state is random if these coefficients Xij are random. The simplest and 
most common random state is to choose Xij as independent and identically dis- 
tributed Gaussian variables [16]. However, the set of complex Wishart matrices 
invariant under every unitary transformation but without any other constraint is 
referred as complex Wishart ensemble or Laguerre unitary ensemble [17) . Its joint 
probability density function (p.d.f.) of N unordered eigenvalue xi,X2, ■ ■ ■ ,xn of 
complex Wishart matrix W is written as 

1 ^ 

(1.4) p^^^-^(xi,...,x^)=- n k-^,f 

l<i<j<N 2=1 

where s > and a = M — N > 0. In this present paper, we extend the scope of the 
index a and s to a > — 1 and Res > 0, respectively. The partition function reads 

N 

(1.5) Zs = ,s^(^+") l[{r {1 + ma + 

i=i 

calculated in the book of Mehta . On the other hand, in case of a random pure 
state 1$), all the eigenvalues ofW = XX'' are not quite same as these of complex 
Wishart matrix due to the additional constraint that trp^^trVF = 1. Thus, the 
eigenvalues of the reduced density matrix pA are distributed according to (|1.4I1 in 
addition to the constraint J^iLi Xi — I. More precisely, 

N N 

(1.6) P'/{x^,...,x^) = -S{J2x.^r) n \^r-^,fY[^T, 

^ i=l l<i<j<N i=l 

where 5{x) denotes the Dirac measure, r = 1 and a ~ M — N . We refer this 
ensemble as fixed trace Laguerre unitary ensemble (FTLUE), following the classic 
book by Mehta [IT] where he referred to fix trace and bounded trace Gaussian 
ensembles as restricted trace ensembles (this class of ensembles has been generalized 
in [1]). We will extend the scope of the index r to Rcr > 0. It follows from (|1.5|) 
that the partition function Z^ equals 

n 7^ ^ NiN+a) H^ll ^1 + j)r(^ + j) 

^ ' ^ T{N{N + a)) 

When M = N, the joint p.d.f. of (|1.6p in |^ is referred to the ensemble of 
random density matrices with respect to the Hilbert-Schmidt metric in the set "Dn 
of all density matrices of size N. It is worthy of stressing that another ensemble of 
random density matrices with respect to the Bures metric is quite distinguished, 
because its features support the claim that without any prior knowledge on a certain 
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density matrix, the optimal way to mimic it is to generate it at random with respect 
to the Bures measure 

The study of the eigenvahies of the reduced density matrix pA = W is crucial for 
understanding and utilizing entanglement. In principle, all information about the 
spectral properties of the subsystem A, including its degree of entanglement, can be 
encoded in the p.d.f. of (|1.6p . For example, one classic measure of entanglement is 
the von Neumann entropy defined by S — — tr Inp^ = —J2i=i^'i l^^ij which is 
a random variable. The average entropy (S) is close to In iV — ^ for large N when 
M > N [53]. Besides, some known results on the FTLUE which are the same in the 
limit as these of the LUE include: the global density [25J[3J, the largest eigenvalue 
distribution [3T], the smallest eigenvalue distribution when M = N 16J. 

In this paper, we focus on the so-called correlation functions of the FTLUE. 
We also consider another closely relevant ensemble: bounded trace LUE (BTLUE) , 
whose joint p.d.f. for the eigenvalues is given by 

N N 

(1.8) P^/(:,,,...^:c^)^—9{r-J2^^) U k-x.fH^, 

S i=i l<i<j<Ar i=l 

where 9{x) denotes the Heaviside step function, i.e., 9{x) = 1 for a; > 0, otherwise 
9{x) = 0. Note that the FTLUE or BTLUE bears the same relationship to the 
LUE that the micro-canonical ensembles to the canonical ensembles in statistical 
mechanics. Section 27 in jTT], Mehta posed the "equivalence of ensembles" problem 
whether all local statistical properties of the eigenvalues between fixed trace and 
unconstrained random matrix ensembles are identical, and further speculated that 
working out the eigenvalues spacing distribution for bounded trace ensembles is 
much more difficult. Although universal local results have been obtained for very 
broad classes of canonical random matrix ensembles [HI [HI [3 [55] j only very few 
results on the local limit behavior of the correlation functions for the restricted 
ensembles (no orthogonal polynomial techniques are available!). Recently, some 
progress has been made for fixed trace Gaussian ensembles [2] , [121 E] ) [E] • Be- 
fore we state our main results, let us first recall that the definitions of correlation 
function and some universal results on the LUE. The n-point correlation function 
R^^^^^'ixi, . . . , a;„) of the LUE is defined as [H] 

(1.9) _R^'^-^^''(a;i, . . . ,a;„) = -— — T / P^^^'" {xi, . . . ,xiq)dxn+i ■ ■ ■ cLxn ■ 

Analogously, the rt-point correlation function of the FTLUE or BTLUE is defined 

as 

(1.10) i?*'''(a;i, . . . ,a;„) = — — T / Pf,''' {xi, . . . , XN)dxn+i ■ ■ ■ dxN, 

[I\ - n)\ Jrn-,. 

where </> denotes S or 9. In particular, when n = 1, r'^^^'^ is called the level density 
or the density of states. A classical result for the LUE says that 

where the symbol l(o,i](2;) denotes the characteristic function of the set (0, 1], and 
ipix) is the Marchenko-Pastur law [TH]. However, for n > 2, the study of a finer 
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asymptotics near a point of the spectrum shows [221 [9]: in the bulk, i.e., u E (0, 1), 
where Ksin{ii,ij) = jg ^i^q so-caUed sine kernel; at the soft edge , 

where 

, , Ai'{u)Ai{v) - Ai'{v)Ai{u) 

U — V 

and the Airy function Ai{x) satisfies the equation Ai"{x) — xAi{x); at the hard 
edge, 

where 

° ' 2(u — 

and Ja{z) denotes the Bessel function of the index a. 

On the other hand, the limit global density of the FTLUE (see ^F' ) and the 
BTLUE is also the Marchenko-Pastur law, i.e., 

1 <p "+° 

A'— foo iV 

where denotes 5 or 6*. In the case of the BTLUE, we will prove the claimed result 
in Sect. 5. Considering universality in the bulk, at the soft and hard edges of the 
spectrum of the restricted trace LUE, we have the same local limit behavior as that 
for the LUE. 

Theorem 1. Let R^^f he the n-point correlation function of eigenvalues of bipartite 
entanglement for a random pure state, defined by il.lO\) . The following asymptotic 
properties hold. 

(i) The bulk of the spectrum: for every u G (0, 1) and G M, 1 < i < n, 

uniformly for ti, . . . ,i„ in compact subsets ofW and for u in a compact subset of 
(0,1). 

(ii) The soft edge of the spectrum: for any f G Cc(K"), the set of all continuous 
functions on R" with compact support, 

f{ti, . . . ,i„) det[KAny{tj,tk)]lk=i dJ't. 
(Hi) The hard edge of the spectrum: 

1 r JV + n f f \ 

A- (16A^^" ' (l6^'-'T6^)=det[i.^^(^»'^.)]".-' 
uniformly for ^i, . . . , t„ in bounded subsets of (0, oo). 
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Theorem 2. Let R^^^ be the n-point correlation function of the bounded trace LUE, 
defined by \l.l(ffj . Let f G Cc(i?"),. the set of all continuous functions on R" with 
compact support, the following asymptotic properties hold, 
(i) The bulk of the spectrum: for every x € (0, 1) , 

lim j^rrr^ f fih,--.,tr,)R'n'^ (x+ x+ d"t 



fitu. ■ ■ , tn) det[i^,in(t„ = l d'^t. 

(ii) The soft edge of the spectrum: 



( Hi ) The hard edge of the spectrum: 

lim [ f(ti,. . . , tn)Rn^ ( . . . , d"t 

fit,,...,t,,)det[KjJt,,t,)]lj^^d^t. 

To the best of our knowledge, this is the first result about the local properties of 
correlation functions for the bounded trace ensembles. Theorems [1] and [5] give an 
afhrmative answer to Mehta's "equivalence of ensembles" problem in the case of 
Laguerre unitary ensemble. In fact, our method can deal with some more general 
ensembles which will be considered in a forthcoming paper. 

The plan of the remaining part of our paper is the following. Sections 2, 3 and 
4 are devoted to the proof of Theorem [TJ Section 5 deals with Theorem [21 In 

Sect. 2, the asymptotic behavior of * in the bulk of the spectrum is given 
based on the rigorous estimates of the correlation function R^^^'"' in the complex 
plane, inspired by [llj . Some of the results in [28' play an important role on our 
proof. In sect. 3, by using the similar method introduced by the authors in |15j . 
the universality at the soft edge of the spectrum is proved. In sect. 4, based 
on a heuristic idea in [2] where universality at zero is considered for fixed trace 
Gauss-type ensembles, the asymptotic behavior at the hard edge is derived. In the 
last section, a "sharp" concentration phenomenon is observed, then local statistical 
properties of the eigenvalues between the fixed and bounded LUEs can be proved 
to be identical in the limit. So we extend the results from Theorem [1] to Theorem 

m 

2. PROOF OF THEOREM [TJ THE BULK OF THE SPECTRUM 

2.1. The relation between i?^^^^*^ and Rf;''. For every i9 e M, through this 
section, we will denote by (.)'' the function 

(2.1) {.f : C\(-oo, l]^C: expi^logz, 

where log denotes the principle branch of the logarithm. The constants C{fj,), Ci(/z), 
C2(/i), ^ifJ-), depending on the parameter /i, may change from one line to another 
line. 
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Let s > 0, and let for every fc > 0, hk{x, s) be a polynomial of degree k with 
positive leading coefficient such that for hk{x,s), 

/>oo 

(2.2) / hk{x,s)hj{x,s)x"e-''^''dx = Sj^k, j, k^O,l,---- 



The generalized Laguerre polynomials {L"{x^ ^))j>o with the positive leading coef- 
ficients ji are defined by 

oo 

(2.3) ^)^' = (1 + '""y"'^ exp(-^), 
^ — ' 1 + sw 

3=0 

and for i, j = 0, 1 . . ., satisfy the relation 

(2.4) e-^'^x^L': {x, s)^ {x, s)dx = ."+^+'+^- ^^^(^ " ^ ^ S,,, . 



Hence 

r(i + i) 

The functions defined by 

(2.6) >f,{x, s) = x"/^e--/^^'%{x, s),i = 0,1,..., 

form an orthogonal sequence of functions in the Hilbert space L^(0,oo). Next let 
us consider the standardized Laguerre polynomials [24' with the positive leading 
coefficient with respect to the weight x"e~^ by the relation 

(2.7) Lf{x) = L^ix, 1), h{x) = Ux, 1), = 1), 

which satisfy 
(2.8) 

i,^(a;,s) ^ s'Lf{xs-^), hi{x,s) = s-'^°'+^^''^h'^{x s'^), g>^{x,s) = s^^/^ ^,(a; g-^). 
The following three recurrence formula for L"{x, s) holds: 

(2.9) jL-^ix, s) = {x ~2j-a + l)sL'^_^{x, s) - s^L^^^i^, s){j + a - 1) 

where j = 2, 3, . . .. The n-point correlation function Rf^^^^'' of the LUE could be 
expressed as 

(2.10) Rn^'^'-'i^u . . . , a:„) = detiKN{x^,x„s))l^^„ 
where 

Af-l 

(2.11) Kn{x, y,s) = ^ <fk{x, s)(pk{y, s). 
Let 

(2.12) KNix,y) = KNix,y,l), 
from (|2.8p . we get 

(2.13) KNix,y,s) = s-'KNixs-\ys-^). 
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The ChrisofFel-Darboux formula for kernels and K^q reads 

(2.14) Kr,{x^y) = ^ ^^^ + ^fNi^)^N-,iy) ~ ^Niy)m-i{x) 

x-y 

and 

(2.15) KN{x,y,s) ^ ^yN{N + a) . 

x~y 

Note that the reproducing kernel Kn{x, y) has the following integral representation 
(Eq.(4.2), [29], or Eq.(3.6), M) 

(2.16) KN{x,y)^ ^^^^ / Si{x + z)S2{y + z) + Si{y + z)S2{x + z)dz 

^ Jo 

where 

(2.17) S,{x) = y]V^ + ViVT^ ^^'^^^\ 

x X 

(2.18) 52(.) = ^/]vT^^ + V^V^^^, 

X X 

and so, for every s > 0, the following relation 



(2.19) Xa.(-,^)^ + r°°5i(-+z)52(^ + z) + 5i(^ + z)^2(- + ^)d2 

holds. The relations p.8|) . ()2.9p . (|2.13p and our extension of the function (.)'' to 
C\(— oo,l] allow us to continue Li{x,s), (pi{x,s) and Kjsi{x,y, s) to this domain 
analytically in the parameter s. The relations (|2.3p - ()2.15p remain valid under these 
continuations. So does (|2.19l) whenever this integral is well-defined. 

Next, for R'^^^''^ and R^^ , we will prove that one can be expressed by the 
other. Let us stress that the similar relation will be frequently used in the proofs 
of Theorems 1 and 2. 

Proposition 3. Let R^^^^'^ and R^^f he the n-point correlation functions, defined 
by il.9\) and il.lO\) respectively, then we have the following integral equation 



(2.20) i?^^^'^(a;i,...,a;„)= / R'^^x,, . . . ,Xr,)j{-)s-^du, 

Jo s 

where j{x) is defined by 

Proof By ^1^1, for any / G L°°(M"), 

f{xi,...,Xn)R^^^'''{xi,...,Xn)dXi...dXn ^ ^' 



ZsiN-ny.Jo J A 

N 

l<i<j<A' 1=1 

Here we make the change of variables: x^ — uy^, I < i < N , and yi belongs to 
the standard N-simplex A = {{yi, ■ ■ ■ ,2/Ar)| J2iLi Vi = 1j 2/i ^ 0}- By (|1.10p . the 
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right-hand side of the above equaUty equals 

zl j 

Zs Jo 



^1 pCO 



■^s Jo JR" ^ U 

= YJ^ /(2;i,...,X„) R^rix,,...,Xn)e-^U^"+^''-'dudxi...dXn. 

Here we have used the foUowing fact that 

(2.22) Rf{''{xi, . . . , = R'/{xir-\ . . . ,x„r-i)r-". 

Thus we prove that 



(2.23) i?^^^.^(2;i,...,^„) = ^ / R':;^ixi,...,Xn)e-^u'' +'"'-'du, 

Jo 

It follows from pTS)) and pTT]) that 

(2.24) = r{N{N + a))s^(^+")Z]. 

This proves this proposition. □ 

Note that the relations Ea. (|2.20p and Ea. (|2.22p also hold when extending the 
parameter s and r to the domain Res, Rer > 0. By Proposition 2, set u = 
{N + a + v)/4:, s = l/4iV, we have 

(2.25) i?r - / i?„ " A^7(iV(iV + a + i;))dv. 

Jo 

We will state the following lemma, which plays a central role in our proof of uni- 
versality in the bulk of the spectrum. 



Lemma 4 

(2.26) / Rrl N-f{N{N + a + v))exp{-iyv)dv = (t)N{y)R'. 
Jo 

where 

(2.27) 0^(y) = e'^(^+"Hl + ^)-(^'+^"). 

Note that the function 0Ar(-) is the characteristic function of A^7 {N{N + a + •)). 
Proo/. Set (iV + a + u)/4 = M, the left-hand side of Ea. ((2?26|) equals 

i?^;''4iV7(47Vu) exp(-iy(4u - iV - a))du 

i?^^"47V (47Vm)^ +^"-1 e"*^" exp(-iy4M)du 



r(Af2 ^ jsTa) J^ 







iy{N+a) ■ 

T{N^ + Nay N' 
Here we have used Proposition O □ 
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Remark 5. By Lemma HI using the inverse Fourier transform, we know that 
(2.28) R'^-^N^iN{N + a)) = ^ </>^(y)i?/ '^"^"'dy. 

2.2. Estimate of R^^^'^ in the complex plane. For the convenience of the 
reader, we wih review some basic results in [28J. The function (Eq.(3.38), ^5] ) 

2 (z- 1)1/2 

ITT 

with principle branches of powers denotes analytic continuation of the standard 
Marcenko-Pastur law in the domain C\[0, 1]. The following two formulas: 

(2.30) g{z) = f log(z - yMy)dy, for z e C\(-^, 1] 



(2.29) ^(z) = - ^ ^^/^ , forz£C\[0,l] 



and 

(2.31) C(z) = -iiT I ^p{y)dy, forz e C\(-oo, 1] 



come from Eq.(3.30) and Eq.(3.40) respectively in |28j. The uniqueness of analytic 
function (Eq.(5.2), gS]) shows that 

(2.32) 2^(z) = 2g{z) -2z-l, for z e C\(-oo, 1]. 

Here I is given by Proposition 3.12 in [2,8] • 

For < < 1/2 and /3 > 0, the sets Sg^p and Sg^p are defined by 

Se.,fi = {zeC\e < Re(z) <l~0, |Imz| < 



Se,p = {zeC\e< Re(z) <l-0, |Imz| < f3}. 



Lemma 6. Let A(i7) = 1 + iH . For every < 6 < 1/2, there exists a positive 
number Hq{9) > 
A > the relation 



number Hq{9) > such that for kernel Kpf{x,y, ^) denoted by i2.11]) . and for 



\ ~ f Uj v 1 

tt{u-v)A{H) 

holds uniformly for all x e [6,1 - 6], < H < Ho{9) and \u\, \v\ < A. 

Proof. For any fixed < 6 < 1/2, taking 6 — 9/A (see Figure 5 in the section 3.8, 
dH]) such that 

[9/2, 1 - 9/2] Pi dUs = 0, [9/2, 1 - 9/2] Q dUs = 0. 

Here dUs and SC/a denote the boundary of the two disks Us and Us around unity 
and zero respectively, depicted in Figure 1 (See Figures 1 and 5, I28j). Hence there 
exists a positive number /3(|) such that this set 5*8, /3(|) is the subset of the set 
^<5 \}Bsy}{5,l — 5), described in Figure 1. Let 

U V 

(2.33) ^ (. + ^)A(i7), = (. + 

It is not difficult to check that for all x E [9,1 — 9] and |u| < A, there exists 
Nq, when N > Nq , Rewi,Refi;2 G [9/2,1 — 9/2]. Furthermore, there exists a 
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Figure 1. As and Bg depend on the parameter 5 and the fixed 
angel v S (27r/3, Stt/S). Details about i/ is described in the section 
3.8 and Proposition 3.19 in [25] . 



positive number Hq{6) such that < Imu'i,Imii;2 < /3(|) for all < < Hq{9) 
if N is sufficiently large. Note that wi,W2 € Ag. From a series of transformations 
yi-^J7h^Th^S'h^i? (See (3.4), (3.13), (3.36), (3.51) and (3.99), [21]), we obtain 

(2.34) 

Kn{wi,W2,^) ^ mKNi-iNwi, 4Nw2) 

~2Nwi -2NW2 (\ 

^\N{4Nw^)'i{\Nw2)^^--^, -(0,l)y-i(4iVu;2)l^(4A^u;i) L 

X 5(z«2)-i5(u;i) J) e("^^(»i)-i"0-3 (47V)-t- Q , 

where (73 denotes the third Pauli matrix. It follows from Eq.(3.99), Eq.(3.56) and 
Theorem 3.32 in [55] that the following relation 

(2.35) S{w2r^S{w,) = I + 0{jj) 

holds uniformly for all wi,W2 G <S'0/2,^i(e/2)- Actually we obtain 

S{w2)-^S (wi) =1 + 0{wi - W2). 
Combining (lOi)) and ([05]) . we find 

-. a/2 a/2 

(2.36) Kn{wi,W2,—) = 2i^lw'^~w ^ [wr"e<(^^)-"^'^^^ -wf "e"^^"'^)-"^^"'^)]. 
Here we have used Eq. (|2.32p . By (|2.33|) . the following asymptotic behavior 

< {w2)-n^ {w,) = -iVTT* / V(2/)dy = \ ^^ A(g)(^-t;)+0(-) 

holds uniformly when x,u,v,H satisfy the assumptions of this lemma. By (|2.36p . 
we conclude the proof of this lemma. 

□ 
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Next, first we will obtain an upper bound about orthogonal polynomial /lAr(z) 
in the complex plane. As a consequence, basing on the integral representation 
Ea. (|2.19l) . we can derive the upper bound estimate of the reproducing kernel i^Ar(-, •, s) 
(See Lemma 8). 

Lemma 7. For every ^ > 0, there exist constants C(^) and f2(u) such that the 
two inequalities 

(2.37) \hNi'iNz)\ <C{^i)N-^n{p)^\z\^+i 

and 

\hM-MNz)\ < C{ii)N-'^ni^i)^-'\zf~'2 
hold for every N and every z satisfying Re 2 > 0, Imz > /i. 
Proof. Theorem 2.4 in shows that for Re z > and Im z > ^, 



hNiiNz) = (47Vz)-te2^' ' ^ 



2Nn 

(2-38) 2zmi^^A)m exp(— A^/' i'is)ds){l + 0(1)) 

where the error term is uniform, and 

(2.39) *(z) = 2z- 1 + 2zi/2(-^_ ;L)1/2^ for z e C\[0, 1]. 
Note that the function g{z) satisfies 

(z- 1)1/2 

(2.40) g'(z) = 2(l-^ > ^ 



z 



1/2 



with the initial condition (See (3.26), (3.28) and (3.35), [28]) 

(2.41) e^ffW =z^ + 0(z^-i). 
A direct computation shows that 

(2.42) g{z) = 2z-l- 2z^/^{z - 1)^/2 -2ln{z^/'^ - {z - 1)^^^) - 2 In 2. 
Note that ? = -2 - 41n2 (Remark 2.3, [H]). It follows from (^311) and ([T^ that 

(2.43) h.i^Nz) . (4A-^)--V2]V^ J/;(;!,)i/.(,i/..(,„,)i/.).. 

For a given analytic branch, {z}!"^)"^ = z, (z — 1) = ((z — X^I'^Y ■ Set w = z^^^, 
which is a univalent analytic function in C\(— oo, 1]. Let x{w) be the root of this 
equation: x + ^ ^ 2w satisfying |x| < 1. Hence 

(2.44) x{w) = w - v^w2 - 1, 

where the function — 1 is a univalent analytic function in C\(— oo, 1] satisfying 
the relation ('v/u)^^^)^ > for w > 1. For every < r < 1, the equation |a;(w)| = r 
defines the ellipse 

^'■''^ (FT7F+(F^^^' 

and |a;(u))| < r denotes the outside of this ellipse E,.. Thus for w G E^, 
(2.46) \x{w)\ ^\w- \/w^ - 1| < r < 1, 
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and so 

(2.47) Re{w - \/ - 1)^ < 
By Eq. (P^ . for w G -E^ we get 

(2.48) min = \i}--v). 

\x{w)\=r I r 

Hence, 

(2.49) }—^ = ^—^\2w-x{w)\ 

(2.50) < 1 + 2\w\ = \w\{2 + ^) < kl(2 + (T^)- 
Further, we have 

\w'-l\ = \w-l\\w + l\ = \lix+^)-l\\^ix+^) + l\ 

(2.51) ^|J_|2|^_l|2|^ + i|2>_l_(i_^)2_ 

Apply (p:iS)) - (P3T|) to (I132D, we find that for every fixed < r < 1, the function 
\hN{ANz)\ could be controlled by 

^ |gAr(.i/^-(2-l)i/2)2|^^ 
2iV^ 2|zl/4(2 _ 1)1/4| 1(^1/2 „ _ 1)1/2)|2 

— 1 |giv(.w)^|gjv 



|;.Ar(47Vz)| = |(4iVz)- ^, .MM,,.. 1/2 lM/2N|2^+a+l 



^ 2iV7r2|u;|l/2|^2_ l|l/4 |a.(y,)|2Ar+a+l 

(2.52) < C,{r)N-^Lirf\wr+^ = Ci{r)N-^ L{rf\z\^+-- . 

Let = n^^- Note that by the convexity of C\Er, e E,. if z e Thus we 
have 

\hN-ii^Nz)\ = |/ijv^i(4(iV- l)zjv)| < Ci(r)(7V- l)-^L(r)^-i|zjvr-^ 

(2.53) < C2(r)iV^'^L(r)^|z|^-3. 

Here we have used the following fact that 

lim ( )" + 2 = e . 

According to the assumption about z of this lemma, it follows from w = |z|^/^e*^^ 
that Raw > 0, Imtu > 0. Note that w'^ = (Rew)^ — (Imw)^ + 2iRewImw. The 
assumption Rez > implies that Re it; > Imw. For every /i > 0, let 



+ 2 I fj, 
2' 

then < r < 1 is a solution of the equation (^ — r)^ ~ 2^. Hence Imz > /i implies 
that 

2(Rew)^ > 2RewImti; = Imw^ = Imz > -( r)^, 

2 r 
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SO Rew > -^{-^ — r), which estabhshes this fact that w £ Er- It follows from 
Eq. (|2.52p and Eq. (|2.53p that for every fi > 0, there exist some constants C(/i) and 
such that 

(2.54) |ftjv(4iVz)| < C{^i)N-^n{fi)^\zf+^, 

(2.55) |/iw_i(4iVz)| < ci^i)N-^n{^i)^-^\z\^-i. 

□ 



Lemma 8. For every /j, > 0, there exist C{ii) and ri(/i) such that 

(2.56) ^\Kn{z^,Z2, ^)| < C{^,M^,f''-'T{2P + 1)1 Imzil^+^l Imzil^+5 

holds for all complex numbers Zi, i = 1, 2 satisfying Rez^ > 0, Iniz^ > /i, and 
N > Here P = N + [^-^] + 1, [2:] denotes the integer part of a real number x. 

Proof By (l^jg)) . one finds that 

1 ~ , 1 , , , JN(N + a) 
-Kn{zi,Z2,—)^ AKn{^Nz^ANz2) = 4^^^^ -^N 

(2.57) 

r' + OO 

Sx{AN{zi + 6i))S'2(4iV(z2 + 6)) + Sx{AN{z2 + 6i))52(4iV(zi + 6l))d6i. 



It follows from Eq.jSj]), Eq. ([2T7ll and Lemma [7] that 

\s,{m{z^ + 9))\ < Vn + VnT^ \miz, + e)\ 

^ ^^ \4N{zi + 6>)|"/^e-^^ "°(^i+^)C(/i)jV-'^»(^)^|^i + 6l|^+5 



|4iV(zi +61)1 

, |4iV(zi + 6l)|"/2e-2Af i^'=(^i+«)C(Ai)7V-'^l)(Ai)^-i|zi + 6*1^-3 

+ mWT9j\ ■ 

Here we have used this fact that 9 > 0, Re(zi + 9) > and Im(zi + 9) > ji. A 
direct calculation tells us that 

|S'i(4iV(zi +9))\< A°'/^-^N-^C{^l)n{^l)'^\zl + 6»|W+^e-2^^'=(^i+^' 
+ A'-l^-^N-^'^^/NTaC{^l)n{^lf-^\z^ + Q\N+^-i^-2NK.(z^+e) 

= ^^n-lN-'c{^,)n{^.r-\-^^^<^^+'^\z^^9f+^^ 

< 7V-iC(Ai)r!(A*)^-ie^2iVRc(.i+e)|^^ ^ 9f+^. 
Note that |zi + 9\^^ < l/fi. Analogously, we have 

(2.58) \S2{'iN{z2 + 9))\ < N-^Ci^^)n{^i)^-^e-^^^''''^''+'^^\z2 + 9\^+^ . 
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Applying the two estimates of Si and ^2 to (|2.57p . we find that 

(2.59) < C{^iM^lf^-^ e-^NRoiz,+e)^^^ _^ 9\^^^+^U0^ ' 

1 

Let {^^} = — [^^]- For every complex number z satisfying Rez > 0, 
Imz > /i, we obtain 

+ 00 Z'+oo 

Jo 

f + oo 
JRcz 

JRcz 

<^2({^}-i) / e-4^''|r;2 + (Iniz)2|^d77 
^0 

<^2({^}-l)(j^^)2P+l / e-4JV,M|^2^^|P^^ 







^2({^}-l) 



< 



+ 00 



2({^}-i)(Iin^)2P+i / e-^A^^^h + lp^dr? 



(2.60) =^2({^}-i)(j^^)2P+ig4iVM( i )2P+i / e-V'^rfy 

4A'/i Jam,, 



Here we have used the inequality: 77^ + 1 < (ry + 1)^ for rj > 0. We also notice that 
(P.17, [TO]) 

1 /■+°° c c™ 

x^e ^dx = €^"{1 + 77 + H ^ r), for to G TV, c> 0. 



r(TO+l)y^ ' 1! 2! to! 

Then we have 

^Nii 74Ar^ AN^i J^Nf. r(2P + 1) 

= e^^-(^)^^+ir(2P 
(2.61) 
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Here we have used the inequahty: AN/j > 1. It fohows from ()2.59|) . p.60p and 
(p:BT|l that 



1 



r(2P + l)(Im (Im Z2)^+^. 



Hence we complete the proof of this lemma. 



□ 



Corollary 9. Let R^^^ '^ he n-point correlation function of the LUE, defined by 
il.9\) . for all real u, ti, i = I, . . . ,n, if N is sufficiently large, the following inequality 



7V« 



N^{u)' 



)A{H) 



< n!C(0)"r!(0)(2^-2)np(2p+ l)«//(2P+l)n(l _ ^)(2P+l)n 

holds in all u e [0,1 - 9] , \t,\ < A and H > Ho{0). Here < 9 < I, A > and 
Hq{9) as in Lemma\^ 

Proof For any u e [9,1- 9], there exists N' such that [ti/Nilj{u)[ < 9/2 for all 
N > N' . Hence we have 



Re((ii 



t,. 



Nij){u) 



)A{H)) > 61/2 > 



and 



Im((M - 

It follows from Lemma [8] that 



< C{9)n{9f^-'^T{2P + l)[H{u - 



NtP{u) 



< C{9)Q.{9f^-^T{2P + l)H'''+'{l - 



2P+1 ; 



\2P+1 



By the definition of determinant and Eq. ()2.10p . we have 



N' 



Nip{u 



NtPiu)' 



NtPiu) 

J_^^« I 



< nl max [—Kn((u+——^—— 



)K{H),u^ 



Nii){uy AN' 



< n!C7(0)"r!(0)(2^-2)np(2p+ l)«//(2P+l)n(l _ p(2P+l). 



We complete the proof of this corollary. 



□ 
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2.3. Proof of Theorem 1 in the bulk of the spectrum. Now, we turn to the 
proof of Theorem [T] in the bulk. 

Proof. From the Stirling's formula r(z) = V^n exp{-z)z^-^/'^{l + 0( j^)) if \z\ 
oo, I arg(z)| < TT, we know that 

(2.62) lim N-f{N{N + a + v)) = ^=6-"^ . 



Particularly taking v — 0, then {N{N + a)) — > l/\/2n as iV — > cx3. A straight- 
forward calculation proves that 



(2.63) lim / (j)Niy)dy ^ V2Tr 

By Eq. (f2?28|) . it suffices to show that 

1 /"^"^ , , LUE 



lim 



sin(7r(ti - tj))^" 

h3 = l 



= V27r det 



7r(ti - tj) 



holds uniformly in u, ti, . . . ,i„ satisfying the assumptions of this theorem. Note 
that 



(PN[-y) ^ (PN[y), Rn ' = Rn 



and 

lim / Re0Ar(?/)(iy = — — . 

It is enough to show that 

^ r-+ca LUE ^. t t 

'si n(7r(t, - tj) )^ 

T^iU - tj) y ^ .^^ 

Making the change of variables: y/N — H, we find 



/•oo 

= lim / Re (/)jv(?/)'iy det 



1 P + OO TTIP 1 + + 



N 



{N^{u)Y Jo 



+ 00 

dH 



Re 



(iViI)i?--^« ((. + ^)A(i/), ...,(.+ ^)AW) A(i/)^ 



9 AT 



^ (/ +/ )=/l+/2 



(ivv(«)) 



Jh 



where H e (O,ifo(^')) fixed and Ho{0) as in Lemma 5. Here we have used the 
following fact that 

(2.64) i?^^^--(xi, . . . , a;„) = Rf^^^'^x.a-', . . .,x„a-')a-- 



BIPARTITE ENTANGLEMENT OF A RANDOM PURE STATES 



17 



holds for any Re cr > and s > 0. On the other hand, let 



J Be(f)N{y)dydet 



fH p+oo 



= {["+ r^)ReMNH)dH det f ^^^^^/^^ ~ = Ji + J2. 

Next we will prove that when N ^ 00, 

(2.65) h-Ji^ 0, /2 -)> 0, J2 -> 0. 

By Lemma 5, for every 0<^<1,A>0 the following relation 

\ / i,j=i 

holds uniformly in all u e [6,1 — 9], \ti\ < A and H £ [0,H]. The uniform continuity 
of the function K{u, H)A{H)" with respect to H means that for every e, there exists 
an H{e) < Ho{0) such that 

(2.66) \K{u, il) A(il)" - K{u, 0) | < e 
holds for all < if < H{e). We choose H = H{e). Note that 

(2.67) \ct>N{y)\ = exp(- ^^^^+"S n(l + ^)), 
and so 

(2.68) iV / \Re(l)N{NH)\dH <2 exp(-^)dj/. 
Hence, for large N the difference Ii — J\ can be controlled by 



|/i- Jil <iV 

+ iV 







f 


Re 


Jo 








r 


Re 


'0 





3 AT I 



'AT 



{NH)K{u,H)A{H)" - K{u,0) 



dH 



r+00 y2 

2 



<e4^ exp(-^)dj/. 



The uniform boundedness of the function K{u, 0) with respect to u € [0,1 — 0] 
implies that as A'' ^ 00, the following relation 



/■+00 f+00 2 

IJ2I < |^(w,0)|iV / [Re(f)N{NH)\dH <[K{u,0)\2 exp(-^)dt/^0 
Jh Jnh 2 
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holds uniformly in all u. It follows from Corollary [9] that 

I/2I < ^V?^!C'(f)"^^(6')^^^~^'"r(2P+ 1)"(1 - i)i^P+^)^ 
tp[u)^ 2 

+00 

I Re [(f>N{NH)A{H)''] \H(^P+^>'dH 

H 

< J!_„!C(0)"f](0)(2JV-2)n(i _ ^)(2P+1)„ 

X r(2P + 1)" e-^^ in(i+^^^)(i + H^)^-H^^P+^y-dH 

Jh 

/'+00 
Jh 

/>+c>o -1 

< e^"^r(2P + 1)" / e-(^^-"(^+i)-i) -^dH 
~ Jh i + H^ 

Jo 1 + i?- 

Here the constant C" depends on 9 and n. On the other hand, we notice that 
r(2P + 1)" = V2TTexp{n{2P + 1/2) ln(2F + 1) - 2P - 1) - o(e"2JVin(2W)-)^ 
thus /2 — 7- as 00. We complete the proof of this theorem in the bulk. □ 

3. PROOF OF THEOREM [TJ THE SOFT EDGE OF THE SPECTRUM 

From Eq. (|2.20l) . we get 

(3.1) Rn'''''^ix,,X2,--- ,x^)^ ,^)i^^(7V„u)7V„d7., 

Jq u u u 

where = N{N + a). Next, We prove a more refined asymptotic result than 
Eq.dHai. 

Lemma 10. Let {b^} be a sequence such that bj\i ^ but Nbfq /^/]nN 00 as 
N — > 00, then we have 

(3.2) / -f{N^u)N^du^ / 7(47VM)4iVdM + 0(e-5(^^")'(i+°(i»), 



where u± ^ ^(l±6Ar). 

Proof. Divide the left hand side of the Eq. (|3.2p into three parts 

/ +/ +/ \'i{N{N + a)u)N{N + a)du. 

Jo Jl-bN Jl+bN j 

First consider JJ*" 7(4A^u)4A^ du. By 7'(4A^u) = 0, we get the maximmn point 

A\ (A^ + a) 

(3.4) Umax = > M_ 

for sufficiently large iV. Note that the function 7(4A^u) is monotonically increasing 
when M g (0,w_], thus 

(3.5) / 7(4iVu)4iVdit<7(Ar„(l-&Ar))Ar„(l-6^r). 



BIPARTITE ENTANGLEMENT OF A RANDOM PURE STATES 



19 



It follows from ([2:2T|) that 
(3.6) 7(A^«(l-fejv))iVa(l"&Jv) = 



(1-6 



N 



liVcp-A'c(l-b]v) 



7v,-^»r(iv„) 

Using Stirling's formula, we get 
(3.7) 

\n[-f{N^{l-bN))N^il-bN)] 

= ln(l - br,) ~ N^{1 - b^) + InN^ - [{N^ - ^) IniVc - + 0(1) 
^ ln(l - b^) + N^bN + ^ InN,, + 0(1) 
= N^hbN - h% + Oib%)] + N,,bN + ^ ln7V„ + 0(1) 
= -^^N + N{N + a) 0{b%) + ^ IniV, + 0(1) 
- t^^[l + 0(M + ^lnA^« + 0(^)] 



-N'bUl + 2^ + o(l) + 0{ — )] = --N'bi,il + o(1)). 

In the last two equalities we used fo^v — >■ and NbN / Vln N — cxj as — > oo. 
Next, we estimate ^^^^{ANu)AN du. Notice that the unique maximum point of 

^{ANu) satisfies 

, _ TV^ + l 

Umax — < 

for sufficiently large N , and v? ^{ANu) is monotonically decreasing when u £ 
oo), thus 

oo poo 

7(4iVu)47Vdu = / u-^{u^-f{mu)4:N)du 

/"OO 

< ul-f{4Nu+)4N / u-2^M = 4iVw+7(4iVw+) =7(Af„(l + 6Ar))iV„(l + 6Ar). 
Similarly, we can get 

1 

2' 

Combing p. 71) and p.8l) . this completes the proof. □ 



(3.8) ln[7(A^„(l + 6jv))^a(l + 6iv)] = -^N^bl{l + o(l)). 



Remark 11. In Lemma [TUl let us take bN — N~'^,k G (0,1). It is a well-known 
fact that the scaling at the soft edge of the spectrum is proportional to N~^^^, 
thus we can choose k > 2/3 and give a very close approximation of correlation 
functions near the radial sharp cutoff point. Then using known results about the 
unconstrained ensembles, we obtain Airy-kernel for the fixed trace ensembles. Such 
arguments can also deal with Bessel-kernel at the hard edge. However, it seems to 
be insufficient for proving universality in the bulk. The main difficulty is that the 
"rate" index k has been rather sharp, in the sense that it cannot be replaced with 
a larger number than 1. 
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Proof. (Theorem [TJ the soft edge) 
For / e Cc(K"), From one finds that 

(3.9) 



((2iV)2/3)«7^„^^ ^' ' " \^ ' (2Ar)2/3' ' (2Ar)2/3 



1 



((2iV)2/3)" 



fih,--- ,tn) 



+ CX3 



1 — 6iv /'l+^iv /"oo 



Rn ^ 



((27V)2/3)n 7j,„ 

/((27V)2/3( w - 1) + wyi, • • • , (2iV)^/^(u - 1) + uy„ ) 

= /i +/2 + /3. 

Here we have used the change of variables: 

(3.10) U = {2Nf/'^{u - 1) + uy^, i = l,...,n. 
Choose — N^'' for fixed 7 G (|, 1). By Lemma [TUl we get 

(^+(4^'---'^+(2i#7l)^"^j/ 7(iV«.)iV.^. 
= 11/11-= (]^0(e-V2(A^''«)^(i+°W)) ^ o(l). 

Similarly, again by Lemma ITOl we have I/3I = o(l). Thus 

(3.11) (EH) =/2+o(l). 
Let 



1 



l+bjx 



((2iV)2/3)" 7jj,„ ■' ^" \^ (27V)2/3' ' (27V)2/3y 
j{Nau)Na du. 
Next we will prove that 

(3.12) lim 1/2 - /2I = 0. 

AT— >oo 
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Since / G Cc(M" ) and 

(3.13) {2Nf'''{u-l)+uy, ""-^ y„i ^ I, . . . ,n, 

we can choose a ball Bn of the radius R in R" centered at zero such that supp{f) C 
Bji and 

{((27V)2/3(u _ 1) + uy,, {2Nf'^{u - 1) + uy^)] 

ivi, • • ■ , yn) e supp{f), i-6jv<M<i + few} c Bji. 

From / e Cc(IR"), given e > 0, there exists some 5{e) > such that \ f{xi, . . . , x„) — 

/(yi, ■■■,yn)\ < e, whenever ||( )-(j/i,...,y„)|| <<5(e),V( ),(2/i,---,2/n) e 

-B_R. On the other hand, there exist Nq independent of (yi, . . . , y„) e Bb, such that 

||((2iV)2/3(„ - 1) + uyu . . . , (2iV)2/3(u - 1) + uy„) - (yi, • • • , y„)|| 

< - + (2iV)2/3) < y^iV-T(i? + (2iV)2/3) < S{e) 

for > Nq. Therefore, V (t/i, . . . , y„) G supp{f) 

(3.14) |/((27V)2/3(u - 1) + uyi, . . . , (2iV)2/3(y _ 1) + uy^) - /(y„, . . . , 2/„)| < e. 
Furthermore, we get 

1^2-^2! < ((2iV)2/3)" 

<e(H-o(l))Cfl,. 

Here is a constant and we have used Lemma[T2]below. Hence the relation (|3.12p 
holds. Combining p.9p and p. lip , we have (|3.9p = ^2 + ''(1) ^or sufficiently large 
A^, more precisely, 

((27V)2/3)»1/^^'''" '*")^"''''"" (1+ (2Ar)2/3'--- ' ^ + (2iv')2/3 ) '^"^ 

= ((^]^ £ /(J/- ■ ■ ' ^Z")^""^ (1 + ■ ■ - 1 + (^) ^"^(1 + oil)) + o(l), 

this proves the anticipated result. □ 

Lemma 12. For any fixed R > 0, let Br he the hall of the radius R in R" centered 
at zero. There exists some constant such that 

(3.15) ((^^ /^^ R'n'^i^ + • • • , 1 + (^)^"y < Cn. 
Proof. Given < (5 < i?, there exists No{R,S) such that when N > No{R,6) and 

{yi, ■ ■ ■ ,yn) e Br, 

n 

(3.16) ^((2iV)2/3(u - 1) + uy,f <{R + Sf 
i=i 

where m e [1 - A^"'', 1 + iV'T], 2/3 < 7 < 1. Let 77 G (0, 1) be a real number and 
(/)(i) be a smooth decreasing function on [0, 00) such that (j){t) — 1 for t e [Q,R + 5) 
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and (f>{t) = for t > {1 + ri){R + S). Set fixi, ■ ■ ■ ,a;„) — (j){\\{xi, ■ ■ ■ ,a;„)||) for 
(xi, • • • , Xn) e M". For N > Nq, we have 

- ((2jV)2/3)n ¥'((2iV)2/3(^ _ 1) + uy,, ■■■ , {2Nf/\u - 1) + uy„) 

R'^'^n^ ••• 1 + — W"?; 

^" (2Ar)2/3' (27V)2/3^°' y 

Multiplying by 'y{Nau)Na and then integrating ()3.17p with respect to u on [1 — 
&7V, 1 + ^at], one obtains 



< 



1 r pl+bN , . 

myW L K(2^)'''(" - 1) + "yi, • • • , (2A^)^/^(. - 1) + 



l-Ojv 



7(JV„„)iV„Rj- .■(l + ^^Jl5^,...,l + ^^|=^)<i„<i", 
Here we make use of Eq. p.lip . This completes the proof of this lemma. □ 



4. PROOF OF THEOREM [TJ THE HARD EDGE OF THE SPECTRUM 

First we prove that the n-point correlation function Rf^^ could be expressed as 
the inverse Laplace transform of Rn^^' , which is slightly different from Eq. (|2.26p . 

Proposition 13. Let Rn^^"^" and R^^f he the n-point correlation functions of 
eigenvalues for the LUE and FTLUE, respectively. Then we have the relation 

B^'Ux, x)- ^'^^"^ r-'^\t-(r^''-^)R^'^^'^ (—X, —x)](r) 

where C^^[h{t)]{x) is the inverse Laplace transform of a function h(t), and Na ~ 
N{N + a). 

Proof For h G L°°(R^), let < h{-) > and < h{-) >s denote that the ensemble 
average is taken in the LUE and the FTLUE, respectively. Consider the integral 

» N N 

I[h]= / h{xi,...,XN)5{r-y2x,)r\xtT\\x,-Xkfd^x. 

•^K" 7^1 tl j<k 

Making the change of variables: Xj = ANyj,j — 1, . . . ,N, we have 

» N N 

I[h] = {4N)^'' / h{ANxi,...,4NxN)S(r~ANy"xAT]xfT]\^3-^kfd^x- 
•^R" 7^1 tl j<k 
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Multiply both sides by e and integrate on r from to oo, we get 
f e-*''/[/i]dr = (47V)^° f h{mxi, . . . ,ANxn) 

Jo ^K" 

N N 

i=l 1=1 j<k 

N N 

= (^)^° / M^^i,---,^^iv)exp(-4ivV:ri) n^TTl^^ -^'^l'^ 



=1 i=l j<fe 



^1/4n( — )'^° <M — •)>• 



Here we have made the change of variables Xj = t~^yj,j = 1,...,N, where 
denotes the principal branch of the power for complex variable t. Using the inverse 
Laplace transform, we have 

I[h] = Z,/,^C-'[i^ f- < h{^-) >](r). 



Notice that 



C-\t-^]{x) = ^e{x),Re{^)>Q. 



The ensemble average < h >s reads 



In particular, taking 



h{xi,...,XN)= f{Xi^,..., 
l<ii<---<i„<N 



then we find that 



/ f{xi,...,Xn)Rfl''{xi,...,Xn)(P3 

f 

-1 27ri J_, 



-ioo+0+ 



X / f{—Xi,...,—Xn)Rn ^" {Xi,. . . ,Xn)d"'X 

_ _r(Na)_j_ r *°°+° rt^-(jv„-n) 

X f{xi, Xn)Rn^^'^ i^Xi, ^X^) d"x 

xRi:^^'^iJ^x„...,±x„)dt. 
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Since i?^*^^'"*" and Rf;^ are continuous, we get 



This completes the proof. □ 

Proof. (Theorem [TJ the hard edge) 

Now we make use of the fa 
can be expanded as fohows: 



Now we make use of the fact that n-point correlation function Rn^^" *™ {xi , . . . , a;„) 



(4.1) 



n n 2N--2 

i?^^^'^(a;i,...,x„) = n-<cxp(-4iV^xO 4S...A}^i' 



2N-2 



= exp(-4iV^.O Y: €LuA 

i=l h,.. .,1,^=0 

where — li + a,i = 1, . . . , n. It follows from Proposition [T51 and Ea. (|4.ip that 



_5 jv + Q 
Rn * {xi , . . . , Xn) 



Nr.-1-rfAT \ 2Ar-2 (N) , AT , " 



(7V + a)^"-i(47V)«^^ f^^^ (AN)^^ 



JVcv-l 



, n 2Af-2 ( N+a v^" NAr^_„_^; i 4 



4 f^""^,..^„=o iAN)^h+n 

^{Na) (AT) i' 

(4.2) 

AT " 2A-2 . n 

(iv„)i:'.+"r(iv„-n-EO 

With rescaling Xi = j^j^, we will deal with the different factors in (|4.2p . First, the 
(^-function term 

i=l 
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as — >■ OO. Since ^ < n{2N — 2 + a), the first term in the sum yields 

(4 4) (1 - y^^^A^c-^-i:':-! 

i—1 

4 ^ t- 1 

(4-5) =l-(iV.-n-^/:-l) ^^ ^^ ^ + ^C^(^i,...,t„) 

i—l 

Here Cjv(ii, ■ ■ ■ ,tn) is uniformly bounded in all ti e (0, • • • , i„ G (0, A] for given 
A > 0. 

Before evaluating the factor containing F-functions, we introduce a lemma about 
ratio of two Gamma functions, due to Tricomi and Erdelyi [3T], see also Copson's 

book [g. 

Lemma 14. Let a > 0, for sufficiently large x the following expansion holds: 

(4.6) -^^^^ = E °^"-^^--;i°"^+^^ i^i°^^)(o). 

x"! (a; — a) s!a;^ 

^ ' s=Q 

Here Bi°'\x) is Norlund's generalized Bernoulli polynomial in x of degree s, and 
Bi"'\o) is also a polynomial in a of degree s. 



Note that there are only finite terms in the sum of the right-hand side of Eq. (|4.6p 
if a is an integer. 

For the convenience, we write 

L,(a) =a{a-l)---{a- s + 1)B(°+i)(0). 

It is a polynomial in a of degree 2s. By LemmalM^ the factor containing F-functions 
reads: 

Since Bq'^\x) = 1, we have LqI^u + J2h) — 1- However, iiJ2h ^ the terms on 
the right-hand side of Eq. (|4.7p 

(4.8) (^^^(" + E0 

approach to a finite non-zero number as iV — >■ cxd, not sub- leading. With the aid of 
multiplication and partial differential operators, we can rewrite 



(4.9) Lsin + h) ti • • • in - LsiY, ^^A) A" ' ' ' 
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Combing (gSl, (|431), (j4J]) and gS]), it follows from (g^l that 

-^—R'^'^iJl Jl^) 

(16iV2)« " 467V2' '16iV2'' 

" (167V2)" (l - (A^o - n - ^ - 1) ^2^^ ^ ^ ^ 

(4.10) 

N-^{N + ayj^^l&' iV2(Ar + a)^^ 16^ 
(167V2)„^^ Z.^^^Z.'='*.^^^Z. s!(iV„)-- '^{'i--.'"}467V2^ 467V2^ 

, 1 it t ^v^iELiM) 

where 



f4iii I ^^b^^^Li^h^ 'v' c^^^ f ^1 yi...f y„ 

^ • ^ s!(iV„)'* (16iV2)« {ii,...,i„}467v2'' 46iV2^ ' 

On the other hand, we know from the expansion (j4.1D for the correlation function 
of the LUE that 



(4.12) 



Af^oo (16iV2)« " 46iV2' 'l6iV2^ 

n 2N-2 

A^ (16^2)n ^'^P I Z^4Ar^ ^{ii,...,(.}46A^2^ 46iV2^ 

1 2A-2 

= lim V c^^' f ^1 y. ( ^" y^ 

A^oo (16A^2)„ ^{ii,...A}46iV2^ 467V2'' 

^ ^ ;i:...,i„=o 

uniformly for ti, . . . , i„ in bounded subsets of (0, oo). Comparing (|4.1ip and (|4.12p . 
the terms where s > in the sum of (j4.11|) vanish. Hence we get 

(4.13) lim I = det[Kj^{U,tj)]l^^,. 

A— s-oo " 
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Thus we have that /2 — >■ as iV — >■ cx). Note that 



Km Ir, — hm „ , , , , , , 

^ ' i—1 s— 



2N — 2 



(16iV2)» ^ {'i,-,'™}^i67V2' ^6^2 

^ ' /i,...,/„=0 



We also notice the following fact: for any sequences {aijili ^-^d if ai 

■ ■ ■ + apf ^ C as N CO, and = O(^), then 

(4.14) |ai6W + a.fef) + • • • + b^^'a^\ < ^ - |6f '| = O ( 1 



1=1 ^ 



where St — ai + a2 + ■ ■ ■ + cii- Thus we have 



|/4| <Cn(27V-2 + a)0(^) = 0(l), 
uniformly for ii, . . . in bounded subsets of (0,cxd). Here we have used the fact 

of Km . 

This completes the proof. □ 



5. PROOF OF THEOREM m 

First we give a representation of correlation functions for the BTLUE in terms 
of these for the FTLUE. 

Proposition 15. Let R^^^ and R^f he the n-point correlation functions for the 
BTLUE and FTLUE respectively, then we have the following relation 

(5.1) i?^(a;i,...,a;„)- f\^u^^-^ ^^Rfi^i"^, . . . du, 

Jq u" u u 

where = N{N + a). 
Proof. It suffices to prove 

(5.2) Ri^^{x^,...,xn)^ r ^u''--^{-rRti^{-x^,...,-x^)du. 
For every w > 0, let 

AT 

(5.3) ^n[u) = {(xi,...,xjv)|^a;j ^u,x.j > 0,j = 1,...,7V} 

be a simplex in , which carries the volume element induced by the standard 
Euclidean metric on R^, denoted by u^-'^daN. For h g L°°(R^), let < h{-) >e 
and < h{-) >s denote the ensemble average taken in the BTLUE and the FTLUE, 
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respectively. From ()1.6p and (jl.Sp . we have 



^ [ ^du [ h{xi,...,XN) Yi \x.i - Xj\'^Y[x°'d(JN 

Jo JA^iu) l<t<j<N t=l 

IT-/" ^ 

— / u^^'^du I h{uXi,...,UXN) Yl - Xjl'^Yl^i^^^ 

^0 iA„(l) l<,<,<Ar i=l 



and 



< h{a-) >s 



^ / ^h{axi,...,axN) Y[ \xi - Xj\'^Y[^id''^N 

JAi^ir) l<i<j<N 1=1 

I r _ ^ 
— '^h{arxi,...,arxN) ~ P II 

"'An(I) l<j<7-<7V 1=1 



N 

i<j<N 

Choose o = - , we get 



(5.4) < h{-) >e= ^ I i-r--' < h(- ■) >s du. 

'^e Jo 

Setting h = 1, we get the ratio of the partition functions and Zg. Substituting 
this ratio, we then obtain 



,u 



(5.5) <h{-)>e^ ^u'"--' <h{--)>5du. 

Jq r 1 

In particular, taking 

(5.6) h{xi,. . . ,xn) ^ ^ /(a;^i, . . . ,a;j,J, 

l<ii<---<i„<Af 

we have 

y (xi , . . . , Xfi^Rj^ {x\^ , . . ^ Xfi) d X 

M —I , f ^ \ nS r / \ m 

-^U " du f{-Xi,...,-Xn)Rn {Xl,...,Xnjd X 

/(Xi, . . . , X„) d"x f ^U^--^ C-TR'/C-X^, -Xr^) du. 



r^" u " u u 



Since i?„ and R^^ are both continuous, we complete the proof. □ 

Next, we notice a "sharp" concentration phenomenon along the radial coordinate 
between correlation functions of the BTLUE and FTLUE. Although its proof is 
simple, the following lemma plays a crucial role in dealing with local statistical 
properties of the eigenvalues between the fixed and bounded ensembles. 

Lemma 16. Let {bpf} be a sequence such that ^ but N'^bpf — > oo as N oo, 
then we have 

(5.7) [' N,,u^"-' du= f iV„u^°-idu + e-^'^"(i+°(i», 
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where U- = 1 — b^- 
Proof. 

N^U^'-'-'du = (1 - bN^" - eA'.ln(l-b„) ^ ^N^{-b^+0(bl)) ^ g-7V^6„(l+o(l))_ 

This completes the proof. □ 

Remark 17. In Lemma [T6| let us take b^ — N^'^,k G (0,2). Since the "rate" 
index k can be chosen larger than 1 while the scaling in the bulk is proportional 
to and at the soft edge of the spectrum is proportional to N~^/^, in principle 
we can prove all local statistical properties of the eigenvalues between the fixed 
and bounded trace ensembles are identical in the limit. Such arguments apply to 
the equivalence of ensembles between the fixed trace and bounded trace ensembles 
with monomial potentials, where we exploit some homogeneity of the monomial 
potentials. 

Before we prove Theorem [21 let us prove the claimed result in Sect. 1: limit 
global density for the BTLUE is also Marchenko-Pastur law. 

Theorem 18. Let r\'^ be the 1-point correlation function of the BTLUE, for 
o-^y f G L°°(K)f]Ciip(M) where the set Ciip{R) denotes all Lipschitz continuous 
functions on R, we have 

(5.8) lim / f{x)-R"{ * {x)dx^ / f{x)ij{x)dx 

where ipi^) ihe Marchenko-Pastur law. 
Proof. By Proposition [T5l 

f{x)-R','^{x)dx^ / N^u^--'fix)--R','^{-)dudx 
N JrJo Nu u 

"^""^ f{ux)^R^;^ {x) du dx 



By Lemma [161 





A^aW °~ /(ua;) — i?i * (x) dudx = Ii + L2. 



7V„it^°-i— " {x)dudx 

Nau^^-^du J j^Rl'^ix)dx 

= ||/||ooe^^'''"(i+°W) ^0 
as N 00. On the other hand, 

N^u^--^{{f{ux)~f{x))+f{x))-R^' ' {X)dudx = l21+L22. 
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Since / S Cup(M.), we have 



for some constant L. Here we have used the fact 

hm / \x\ — * {x)dx— I \x\'ih{x)dx. 



Again by Lemma [TCI 



hm I22 = / J{x)%l){x) dx. 

N—>oo ' 



This completes the proof. □ 

Now we turn to the proof of Theorem [51 

Proof. The proof is very similar to that of the soft edge of the spectrum in Theorem 
[U we only point out some different places in the bulk case. 

In Lemma [TBI choose bjy — iV^", k G (1, 2). The change of variables correspond- 
ing to ([3T0[) reads: 

(5.9) ti = {u — l)N xip(x) + uyi, i — 1, . . . ,n 

where fixed x G (0, 1). The condition that — N^'^, k e (1, 2) ensures {l — u)N < 
jy-K+i Q g^g ^ — y -fQj. ^ g Qj^ ^j^g other hand, by Theorem [1] the 

following fact similar to Lemma [T^ is obvious: for any fixed R> 0, 

(5.10) JWTT^ I ^"'^(^ + 1^77H'--- /vTh)^"2^^^«- 
(Nij:[x)Y J Br N-)p[x) Ni}{x) 

Here Bb, is the ball of the radius R in M" centered at zero, and is a constant. 

Using Proposition [15] and Theorem [l] we complete the proof after a similar 
procedure. □ 
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